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Abstract 

We study the pp — > ppn^ reaction near threshold based on heavy-baryon 
chiral perturbation theory. We show that the two-pion-exchange diagrams 
give much larger contribution than the one-pion-exchange diagram which is 
of lower chiral order in Weinberg's counting scheme. We also discuss the 
relation of our results to the momentum counting scheme. 



The near-threshold pp^ppn^ reaction has been attracting much theoreti- 
cal attention, ever since experimental data of extremely high quality became 
available [1]. The heavy-baryon chiral perturbation theory (HBxPT) offers 
a possible systematic approach to the investigation of this reaction. A mo- 
tivation of this study may be stated in reference to the generic NN^NNn 
processes near threshold. Although HBxPT presupposes the small size of 
its expansion parameter Q/A^, the pion-pro duct ion reactions involve some- 
what large energy- and three-momentum transfers even at threshold ^ ~ 
m^mAr). Therefore, the application of HB^PT to the NN^NNn reac- 
tions may involve some delicate aspects, but this also means that these pro- 
cesses may serve as a good test case for probing the limit of applicability 
of HBxPT. Apart from this general issue, a specific aspect of the pp^ppir^ 
reaction makes its study particularly interesting. For most isospin channels, 
the NN^NNtt amplitude near threshold is dominated by the pion rescat- 
tering diagram where the ttN scattering vertex is given by the Weinberg- 
Tomozawa term, which represents the lowest order contribution. However, 
a quantitatively reliable description of the NN-^NNir reactions obviously 
requires detailed examinations of the corrections to this dominant amplitude. 
Meanwhile, since the Weinberg- Tomozawa vertex does not contribute to the 
pion-nucleon rescattering diagram for pp^ppir^, this reaction is particularly 
sensitive to higher order contributions and hence its study is expected to 
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provide valuable information to guide us in formulating a quantitative de- 
scription of all the NN^NNn reactions. 

At threshold the strong (and Coulomb) pp initial and final state interac- 
tions have to be considered in these reactions. In our DWBA evaluation the 
nuclear transition operators are derived using HBxPT, whereas the initial 
and final state nuclear wave functions are calculated in the standard nuclear 
physics approach (SNPA). A serious problem encountered in performing such 
a hybrid xPT (for short called EFT* below) calculation of pion production 
is that the calculation involves uncomfortably high momentum components 
which are present in the nuclear wave functions. The occurrence of these 
high-momentum components goes against the tenet of xPT, which presup- 
poses the existence of a momentum cutoff scale, A^ ~ 1 GeV/c. The high 
momenta arise from two sources. The first source is the large momentum 
components contained in the distorted initial and final wave functions gener- 
ated by a so-called high- precision phenomenological NN-potential, Vn^; see 
e.g. Ref. [2](b). The second source of the high momentum behavior origi- 
nates from higher powers of momentum terms which appear in the transition 
amplitudes generated by higher xPT diagrams [3,4]. Below we will focus our 
discussion on these transition amplitudes. 

In order to eliminate from the A^A^ wave functions the high-momentum 
components that lie above the original cutoff scale of xPT, a suitably pa- 
rameterized cutoff factor is introduced (we make certain the observables are 
independent of this cut-off). This is admittedly an operational remedy, the 
foundation of which needs to be examined from a formal point of view. It 
is also informative and of practical value to examine the use of the "low- 
momentum regime NN potential", Viow-k [5^7]. Viow-k is derived from Viviv 
by integrating out the high- momentum components contained in Vnn- Since 
Viow-k by construction is free from high-momentum components, its use in 
an EFT* calculations for pion production should alleviate the "high momen- 
tum problem" that plagued the past DWBA calculations. From a purist's 
point of view this may not be a totally satisfactory approach but we believe 
that this "pragmatic" method still has its merits. We remark that, as is well 
known, Viow-kS generated from any realistic phenomenological potential lead 
to practically equivalent half-off-shell NN K-matrices and hence the same NN 
wave function. 

We derived the TPE transition amplitude operators [3] using Weinbergs 
chiral counting scheme with expansion parameter e ^ mjr/mjv — 0.15. 
We isolated the high-momentum components of these amplitudes using an 
asymptotic expansion, see Refs. [8, 9]. Hanhart and Kaiser [4] used the 
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momentum counting scfieme (MCS) [10, 11] to evaluate TPE diagrams for 
tfie reaction NN^NNn. The MCS lias tlie expansion parameter e ~ 
(m^/mjv)''^/^ ^ 0.39. Unlike Weinberg's chiral counting a subtlety in MCS is 
that loop diagrams of a given order z/ in e not only contains a contribution of 
order u (the "leading" part) but, in principle, can also involve contributions 
of higher order in e ("sub-leading" parts). Hanhart and Kaiser evaluated 
the "leading" part of the lowest MCS-order TPE diagrams and showed that 
the "leading" parts of the two-pion exchange diagrams, when summed up, 
cancel among themselves, see also Ref. [12]. Wc [14] identified the "leading" 
part of our TPE operators [3] and confirmed this cancellation [8,9]. As will 
be discussed below, we have found however that the remainder, or the "sub- 
leading" parts, of the TPE amplitudes can be at least as large as the one-pion 
rescattering amplitude [8,9]. We consider it important to re-examine the be- 
havior of MCSs "sub-leading" parts of these TPE diagrams in order to see 
whether they can still be as large as indicated by the phenomenological SNPA 
success of the Lee-Riska heavy-meson exchange mechanism. In a forthcom- 
ing publication [13] we will also consider other chiral correction amplitudes 
including the contributions from counter-terms needed to regulate the UV 
behavior of the TPE loop diagrams. 



Figure 1: The two-pion-exchange loop-diagrams discussed in the text. 



The one-pion loop pp ppn^ transition operators were evaluated analyt- 
ically by Rcf. [3] using HB^PT. When these operators are sandwiched with 
phenomenological determined distorted initial and final NN wave functions. 
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we find that the momentum integrals convergence very slowly [8,14]. This 
slow convergence can easily be understood when we adopt the threshold fixed 
kinematics approximation {FKA) to evaluate the amplitudes. We impose the 
FKA on the analytic expressions for the transition operators of the different 
TPE diagrams are given in Ref. [3] and make an asymptotic expansion in the 
two-nucleon momentum transfer k, i.e. \k\ = \p — p'\ —>■ oo. The transition 
operator matrix T of the TPE diagrams is of the form 

T = (^^^ {E.k)t{p,p',x) (1) 

where x = p ■ p' . The generic asymptotic behavior for t{p,p', x) is [14]: 

t{p,p', x) ti {9Al{^fl)f \k\ + t, ln[AV|^n + t3 + x) , (2) 

where t^ is asymptotically /c- independent, and 6t{p,p', x) is 0{k~^). For each 
of Types I ~ IV, analytic expressions for tj's {i = 1, 2, 3) can be extracted [14] 
from the amplitudes T given in Ref. [3]. The first term with ti in eq.([2]) is the 
"leading" term in MCS discussed by HK [4], whereas the remaining terms, 
which we refer to as the "sub-leading" terms, were not considered by HK. 

Table 1: For the four types of TPE diagrams, K= I, II, III and IV, the 
second row gives the value of ti defined in eq.(l2]), and the third row gives 
the ratio Rk = T^/Tjiesc, where Tx is the plane- wave matrix element of T 
in eq.([T]) for Type K, and T^esc is the lowest chiral order one-pion-exchange 
rescattering (Resc) amplitude. The last row gives = Tj^/Tnesc, where 
Tj^ is the plane- wave matrix element of T in eq.([l]) with the ti term in 
eq.(l2]) subtracted. 



Type of diagrams : K = 


I 


II 


III 


IV 


Sum 


{tl)K OC 





-1 


-1/2 


3/2 





Rk 


-.70 


-6.54 


-6.60 


9.19 


-4.65 


R-K 


-.70 


-0.82 


-3.73 


0.60 


-4.65 



In Table 1, the second row shows that the "leading" parts in the MCS 
OC ti of the TPE diagrams I ~ IV sum to zero. This confirms the finding 
of Hanhart and Kaiser [4]. In the third row marked Rk, and using FKA, 
we give the values of the ratio of the transition amplitude in Ref. [3] to the 
rescattering amplitude in the plane wave approximation for the four TPE 
amplitudes. We note that the sum of the TPE amplitudes are non-zero. The 
magnitude of the sum reflects the size of the "sub-leading" parts of the TPE 
amplitudes. In the fourth row we have removed the "leading" term, ti from 
t{p,p', x) in Eq.([T]) and then evaluated the similar ratio R^. We observe that 
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tfie sum of tlie four modified TPE ampiitudes (II, III and IV) is larger than 
the rescattering amplitude [8] . This can also be inferred from the Rk row of 
the table where the ratios of the amplitudes 11:111:1V are about 

-2:-2:3, whereas the ratios of the corresponding amplitudes "leading" terms 
are -2:-l:3. However, the evaluations of the different diagrams appear to 
confirm that the magnitudes of the different diagrams follow the momentum 
counting rule as indicated in Table 11 of Ref. [15]. 

Next we investigate the behavior of the TPE diagrams as we go beyond 
the plane-wave approximation by using distorted waves (DW) for the initial- 
and final-state NN wave functions. For formal consistency we should use the 
NN potential recently derived from HBxPT [16,17], but as discussed earlier 
we adopt here EFT*. As argued earlier, in order to stay close to the spirit 
of HBxPT, we introduce a Gaussian momentum regulator, exp(— p^/A|.), in 
the initial and final distorted wave integrals, suppressing thereby the high 
momentum components of the phenomenological NN potentials, e.g. the 
Bonn [18] or Nijmegen potential [19]. We require that Kg be larger than 
the characteristic momentum scale of the NN NNn reactions, \p\ ~ 
^/m^^mN — 360 MeV/c, but it should not exceed the chiral scale ~ 1 
GeV/c. Part of the work to implement the idea of utilizing a low- momentum 
regime NN potential has been published [20] where use was made of Viow-k- 
One issue regarding the use of Stony Brook's Viow-k in the present context is 
that it is derived with a rather low value of the cutoff parameter, Pmax = 2 
fm"-*^, which is very close to the threshold momentum for the NN^NNtt 
reactions. We therefore have extended the Viow-k potential to cases where 
the momentum cut-off Pmax is allowed to be larger than the original Stony 
Brook value up to 5 fm~^ (this value corresponds to the chiral scale, A^ ~ 1 
GeV/c^). This extended Viow-^k has been used in our recent study of the 
two-pion-exchange (TPE) amplitudes for the pp^ppn^ reaction [9] . 

To evaluate the pp ppir^ reaction at threshold using HB^PT, the 
impulse approximation (LA.) and the one-pion-exchange (Resc) diagrams 
are the lowest order amplitudes (diagrams) according to Weinberg count- 
ing. However, the typical momentum for this reaction at threshold is p ~ 
^J^rrQm^^ which implies that for this reaction we have to take Weinberg chi- 
ral counting with a grain of salt. It was shown early on that in HB^PT 
the LA. and Resc amplitudes interfere destructively resulting in a very small 
cross sections [2, 10]. We therefore expect sizeable contributions to this re- 
action from the TPE diagrams. In Table 2 we show examples, taken from 
ReL [9], of DWBA evaluations for a typical energy Tiah — 281 MeV. Since 
the ti terms add to zero and to improve the numerical convergence, we drop 
the t\ terms in our calculations as was done in the fourth row of Table 1. 
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Thus, in Eq.([T]), we use t*{p,p', x) instead of t{p,p', x), where t*{p,p', x) is ob- 
tained from t{p,p', x) by suppressing the ti term in Eq.(l2]). The partial- wave 
projected form of t*{p,p', x) in a DWBA calculation is written as: 

— j J p'^dp p' Mp' J dx V^iso (p) f{p, p', x) {p - px)^3p^ (p) (3) 
In Table 2 we show the values of the J amplitudes for each TPE diagram. 

Table 2: The values of J, Eq.Q, corresponding to the TPE diagrams of 
Types I ~ IV, evaluated in a DWBA calculation for Tiab = 281 MeV. The 
column labeled "Sum" gives the combined contributions of Types I ~ IV, 
and the last column gives the value of J for Ivr-Resc. For the Nijm93 poten- 
tial case, the results for three different choices of Ag (in MeV/c) are shown. 
For the case with Viow-k, CD-4 (CD-5) represents Viow-k generated from 
the CD-Bonn potential with a momentum cut-off Aiow-k = 4fm~^ (5fm~^). 
The last row gives the results obtained in plane-wave approximation. 





I 


II 


III 


IV 


Sum 


Ivr-Resc 


14ijm : Ag = 600 


-0.12 


-0.12 


-0.57 


0.07 


-0.74 


0.20 


l^Nijm : Ag = 700 


-0.12 


-0.11 


-0.57 


0.06 


-0.74 


0.21 


\/Nijm : Ag = 800 


-0.12 


-0.11 


-0.55 


0.04 


-0.74 


0.22 


^o^-fc (CD-A) 


-0.12 


-0.09 


-0.46 


0.03 


-0.65 


0.23 


Viou,-k (CD -5) 


-0.09 


-0.06 


-0.30 


-0.01 


-0.46 


0.22 


Plane — waves 


-0.06 


-0.07 


-0.30 


0.05 


-0.37 


0.080 



Table 2 shows that the DW amplitudes from the TPE diagrams are only 
roughly of the order of the one-pion rescattering amplitude tabulated in the 
last column. Evidently, when we compare to the plane wave amplitudes, the 
DWBA treatment does affect the relative magnitudes of the various diagrams 
differently. In the last row we show the plane wave amplitudes and we note 
that the "sub-leading" part of diagram III is a factor three or more larger than 
the other amplitudes. The MCS indicates that the "sub-leading" parts of the 
TPE amplitudes should be of the same order as the one-pion rescattering 
amplitude. Furthermore, the expansion parameter in MCS, e ~ 0.4, ideally 
speaking should be the ratio of the amplitudes from the different orders in 
the MCS. When comparing the results in Table 1 and Table 2, we find that 
the "leading" part of diagrams II and IV are almost an order of magnitude 
larger than their "sub-leading" parts. Clearly, as seen in Table 2, we find, 
as expected in the MCS, that the amplitudes from the "sub-leading" parts 
of diagrams II, IV and from the one-pion rescattering diagram are about the 



6 



F. Myhrer 



Two-pion-exchange effects in pp — > ppvr^ reaction 



same magnitude. However, diagram III (the "cross-box" diagram) appears 
not follow the expected MCS behavior. The plane wave amplitude for the 
"sub-leading" part of diagram III is less than 50% than the "leading" part of 
diagram III. Moreover, the amplitude of the "sub-leading" part of diagram 
III is more than a factor larger than what is expected in the MCS. One 
explanation could be that we evaluate diagram III using HBxPT's heavy 
nucleon propagators and not the nucleon propagator which is advocated for 
the MCS [15]. This issue will be resolved in the near future. A final note, 
the lowest- and next-chiral-order (diagram VII in Ref. [3]) one-pion-exchange 
diagrams were found to be same order of magnitude as expected in MCS. 

We have demonstrated [9] that, as expected, the two-pion-exchange loop 
diagrams give very large contributions to the pp ppm^ reaction at thresh- 
old, and that these diagrams will give important contributions to other 
NN — > NNt[ reaction channels. 
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